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Abstract 



^^ ' It is shown that the motion of a muhielectron atom in an external gravita- 

f^ , tional field in a good approximation is described by system of the Mathisson — 

Papapetrou equations, if we put as a classical angular momentum of the atom 
the expectation value of the operator of the full angular momentum of the sys- 
tem, which includes spins of the nucleus and electrons, and orbital momcntums 
^^ , of the electrons in the atom. 
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1 Introduction 

The interest to the problem of the atom motion in an external gravitational fields is 
caused by the following two reasons. 

First, the question about the possibility of using the atoms as the detectors of 
the gravitational fields represents the significant interest for modern astrophysics 
[Q-Q]. Second, because the atom is an extended quantum system, with the inner 
structure (the motion of extended bodies in general relativity is widely discussed 
in the literature), the problem of atom motion has considerable theoretical interest 
too. 

The analysis of this situation shows, that the influence of an external gravita- 
tional field on the atomic spectra can reach the experimentally measurable values 
only in extremely strong gravitational fields (with characteristic radii of curvature 
<10^'^ cm) |5|,|6[g|, or in the case of ultrarelativistic motion of atom as a whole [p|,|lO|]. 
Apparently, in both cases the motion of an atom cannot be considered a priori as 
a geodesic one. So the finding of the equations of the atom motion in an arbitrary 
curved space-time is the important problem. For the case of the one-electron atom 
this problem was solved in the paper Q. Following this work we shall generalize 
received results for the case of multielectron atom. 

For essential simplification of the subsequent calculations, we shall use two as- 
sumptions. First, because nucleus mass // is much greater than electron mass m 
(/z^m), the atom motion can be considered in a good approximation as a motion 
of a classical point particle with the law of motion 

x^=r(r) (1) 

where r is the proper time: dr^ = —c~'^ds^. Moreover, the assumed identity of atom 
and nucleus masses permits the consideration of the atomic motion as a nucleus 
one. Second, we shall describe the motion of electrons in the atom by quantum 
mechanics. Moreover, because the relative motion of electrons in the atom is the 
essential non-relativistic (ii <C c) , so in the special comoving reference frame the atom 
can be described in frameworks of non-relativistic (or quasi-relativistic) quantum 
mechanics p]. Subsequently we shall assume that the atom is in a quasistationary 
quantum state. 

Thus, one can obtain the equation for ^*(t) from the following suggested model. 
The nucleus, in the classical mass-point approximation with intrinsic angular mo- 
mentum (spin), is moving in an external gravitational field and interacting with the 
electromagnetic field of the electrons. The electrons, which for theirs part is moving 
in external gravitational field and in the electromagnetic field of the nonexcited and 
electrons, are considered in a quasistationary bound quantum state. 

We notice, that this assumptions are usual. It allows us instead of the relativistic 
problem of many bodies, which even in fiat space has a lot of difficultiesF], consider 
the motion of electrons in the given electromagnetic field, created by a nucleus, and 
also in an external gravitational field. 

As it is shown in works [0,|^, the generalized Ehrenfest theorem was proved 
for a Dirac particle in quasi-classical approximation. According to this theorem. 



^Quantum mechanical problem of two bodies in a weak gravitational field was discussed in 
work 0. 



the expectations values of the position and spin operators of a electron fulfill the 



Mathisson — Papapetrou equations ||12|,|13|, which are describing a motion of a classi- 



cal particle with a spin in an external gravitational field. Applying this results (valid 
for any Dirac particles) to the nucleus is represented justified. Therefore, to describe 
the nuclear motion in our model, we will use the system of Mathisson — Papapetrou 
equations. For this purpose, in the first equation of this system we have to add a 
Lorentz force, which describes in first approximation the interaction of the nucleus 
with the quantum state averaging electromagnetic field of its electrons. Here, we 
neglect the terms, which contain the second derivative of 4-velocity and describe 
interaction between the nuclear magnetic moment and the field of electrons [p!^ ], 
because the numerical estimations show, that these terms are smaller than all other 
terms of the equation. In result we receive the following system of the equations: 

Dv} 1 : , j.,_^ Ze 



/X-— = -R' ,fci7x%'='P« V, + — F^%„ (2) 

Dt 2c c 

^)-f^ = -^'^n-^, 6)u"5„ = 0, (3) 

Dt c^ Dt 

where Ze, u* = -^ and 5^ denote the charge, the 4-velocity and the classical 
spin of the nucleus respectively, here under the classical spin we understand the 
expectation value of the nuclear spin operator; e^^Pi is the Levi-Civita pceudotensor 
with £^^^^ = {—detgij)~^'^; R^ jki is the Riemann curvature tensor; F*j is the 
tensor of the electromagnetic field of electrons; c is the velocity of light. Here 
and everywhere the latin indexes run from 1 to 4, and greek ones from 1 to 3. 
The signature of the space-time metric is +2. We shall notice, that equation (^b) 
represents the additional Pirani condition [|l5| 0. Due to this condition it is possible 
to enter a 4- vector of a spin Sk and constant fi, which we shall identify hereinafter 
with nuclear mass. 

In the equation (g) of nuclear motion it is unknown only one quantity - the tensor 
of the average electromagnetic field F^^ , created by nuclear electrons. Thus, the 
problem about finding of the equations of motion of atom in an external gravitational 
field is reduced to calculation of the electromagnetic field of nuclear electrons. The 
most simple this problem can be decided in the special, comoving atom reference 
frame. 



2 The Reference Frame of the Single Observer. 

As a comoving reference frame we shall choose a reference frame of the single observer 
|17| . From our point of view, this reference frame is the most convenient for the 
description of objects of the small sizes. 

The reference frame of the single observer is defined by a motion of a single mass 
point ("the single observer"). The world line of this mass point is named basis line. 
Along this line we establish a comoving tirade HIj^Jt), determined by the condition 

/i/.N = -M* accurate to three-dimensional rotations. If we will use the world line of 
the nucleus (||) as a basis line, we will obtain a comoving reference frame for the 
atom. 



^By force of assumptions made above, the Pirani condition coincides with the Tulczyjew-Dixon 
condition @0. 



Three-dimensional physical space in the reference frame of the single observer 
we shall define as a geodesic spacelike hypersurface /, which is orthogonal to the 
basis line in the given moment of proper time r. At each point P, laying on a 
hypersurface P G /(t), we shall put in conformity three scalars X^°^' = aphf^ k^, 
which will hereinafter play a role of the three-dimensional coordinates. Here ap 
is the value of the conical parameter a at the point P, defined along a spacelike 
geodesic in /; A;* is the tangent unit vector to that geodesic at the point a = 0. 

For nonrotating frame (i.e. when the tirade /i/^-, (r) is set along a basic line with 

the help of Fermi — Walker transport) the quantities {X^'^\ct) correspond to the 
Fermi normal coordinates, defined by Synge [0, |l8|. In general case of rotating 
reference frame, the coordinates defined in the same way 

X° = X("), X~^ = CT, (4) 

we shall name as rotating Fermi coordinates |jl9|. In these coordinates the metric 
tensor can be presented as: 



where 



H^m = -^«HM(/3)H^^"^^^'^^ + Oip'), (6) 

e(Q,)(4) = -e(a)(<j)(r)-'^ '^ +6'(a), (7) 

^(4)(4) = -2(^TV(,)XW + 0), (8) 

Oia) = ^i?(,)(^)(.)(4)^(^)^('^) + 0{p'), (9) 

ii(4)M(4)H - ^%)(^)(4)HM^(")) X(^)XM + 0(/). (10) 



Here, we used the following notations: r}{i){j) = diag{l, 1,1,-1) is the Minkowski 
tensor; 

W^^) = ^-)^ (11) 

Dt 

and 

..(") = ^e(-)('')M/.(.),-^ (12) 

are the acceleration and the angular velocity of the reference frame respectively; 
^mwy'r) ig the three-dimensional Levi-Civita symbol; w-'^' (m){n){k) ^i^d 

R'^'' {m)(n)(fe)(0 = K h\m)hln)hik)h\l)R'' bcd;p (13) 

are the tirade components of the Riemann curvature tensor and its covariant deriva- 
tive respectively, determined along the basic line; p = 4/X(„-)X("). 

Because of the small sizes of atom, the condition £(i)(j) <C 1 is allowable even 
for very strong (from the macroscopic point of view) gravitational fields, therefore 
for further calculations we shall be limited by consideration only linear terms in the 



Taking into account that u* = (0, 0,0, c) and gij{0) = V{i){j)^ ^^ receive the 
fohowing expressions for the system of Mathisson — Papapetrou equations in the 
rotating Fermi coordinates (^): 



d_ 



\4){u)(a){f^) 



S{t) + ZeF(^a){A), 



S{a) = e(a){K){T)^^''' S^'^' , '5(4) = 0, 



where all quantities are taken along the nucleus world line {X^°'' = 0). 



(14) 
(15) 



3 Calculation of Electromagnetic Field of Electrons In- 
side the Atom. 



To find how the electromagnetic field of electrons influences on the atom motion, 
it is necessary, first, solve the Maxwell equations in the Fermi coordinates for the 
electromagnetic field, created by electrons in the arbitrary point C^"\ Second, cal- 
culate the expectation value of that field, using the assumption, that the atom is in 
a quasistationary quantum state. 

Solving Maxwell equations 



det(5^,)5""=r(A,fc-A 



k,h 



0, 



(16) 



where A^ = (A(cf),^4 = —^) is the electromagnetic potentials, we will interest only 
quasistationary solutions, i.e. we shall search the field of charges e, resting in the 
points X)^ (^ = 1, 2, ..., Z). Avoiding in this case the consideration of the motion 
of electrons we neglect the effects of interaction of the magnetic field of electrons 
with the nucleus magnetic momentum. These effects are small enough because of 
the small velocity of electrons within the atom (w <C c) . 

Writing down the Maxwell equations in the linear approximation in the £(«)(,) 
and neglecting the terms, containing the derivative on time of the potentials, we will 
have the following expression for electromagnetic potentials of the system of the rest 
charges in some arbitrary point Q"': 



.(c<"')^E£ 



1 



1 + ^^(/3)<^ + ^^(4)(.)(4)(.)C('^) {2&^ - X(^)) + 



mA'') )rW )r(.^) 



+7r^Rmr)Mi,)0'>C''>x)i>x-:' + 0{p-X) 



6R 



(17) 



1 



+ 2^^"^ (4){^)(r) 



(t)H-^ o; ^n (iy)(4)(T)<, c, + 



<^HXW + IflW (,)(,)(4)XHXW + 0{p^: 



Xa(„\X 



(a) 



A{a)^A 



(18) 



where ^"^ = (((") - X]^^), Ra = ^JRTRAia), PA = 

Thus, we also assume, that inside of our system there is no substance or other 
strong fields, resulting to the additional curvature of space time, i.e. the Ricci tensor 



inside of atom is equal to zero {Rij = 0). Otherwise, the presence of the additional 
interaction would generate the effects, which exceed in magnitude the influence of 
the gravitational field on the atom motion. 

Then, in comoving rotating Fermi coordinates the components of the tensor of 
the average electromagnetic field will be equal to: 



-^{q)(4) = (-^(a){4)/ , (19) 



where 



Fi 



H{4) - QHa) _ ^ c\ dr 



V 5C(-) ;^c.)=o ^V 9r )^,.,^, 



z 



E ( ^ + 1: {Wmm.)X'A' + ^w,^ . (20) 



A=l 



p\ PA \r^^}y''>^''}y->-A ' 2c2 



Here, ^^(q.) is the position operator of ^4— th electron; pA = J Xj^i^^-^X^ and the 

operators (p and A are obtained from (|l7| ) and (^) by changing X^ — > X]^ , 
PA — > PA and so on. 

In order to determine the expectation value (19) of the electromagnetic field 
tensor we will use the quasirelativistic two-component representation of the covariant 
Dirac equation. The detailed explanation of this representation, and the consistent 
formulation of the quantum mechanics in the arbitrary reference frames and in the 
external gravitational fields can be found in the book ^. 

According to this approach [P, the covariant Dirac equation can be interpreted 
in general case as a special coordinate representation of the traditional quantum- 
mechanical equation of motion in Hilbert space. Thus, the problem about solving 
of the Dirac equation in the curved space-time can be lead to the similar problem 
of the traditional quantum mechanics in the fiat space - time, and the properties 
of the reference frame and the curved space-time can be taken into account in the 
Hamiltonian as the additional potentials. 

In particular, in quasirelativistic approximation our system can be described by 
a two-component Pauli-type equation 

ih^ = HiP (21) 

OT 

(with 

fiP+i:d^X = l, (22) 

/ 

and d^X = dX''^^ dX^"^^ dX^^^) in the reference frame of single observer. The complete 
expression for the two-component Hamiltonian H (for one-electron atom) can be 
found in ^. The analysis of the complete operator H shows, that the infiuence of 
external gravitational field will always be more significant than effects caused by the 
non-Coulomb terms of the nuclear potential^. 

Therefore, in a good approximation the Hamiltonian of the system of Z electrons 
is given by: 
H = Ho + Hi, (23) 

^The electromagnetic potentials of the nucleus can be obtained from ([17|) and ([IS) by the fol- 
lowing substitution: e -^ Ze, (,% ^^ ^a ' ^^'^ Ca ~^ 0. 
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where 

is the Hamiltonian of the multielectron atom in the flat space-time; 



^(^) (LAm) + SAm) + "i^rm^f 1 (25) 



A 



is the Hamiltonian, describing the interaction of the electrons with the external 
gravitational field. Here we introduced the following notations: X^ = Xj^ , 

PA{a) = jd/dXX , LA(a) = e(a)(/3)(r)-^)f ^i ^nd SA(a) = ^cr{a) are the opera- 
tors of the position, momentum, angular momentum and spin for the ^-th electron; 

A = 1,2, ...,Z; a(^a) are the standard Pauli matrices; 9A{a) = 3P{a){ti){u){4)^A^Ai 

PAB = y{XA{a) - XB(a))iXA " ^B )• 



From (|20| ) by using the symmetry of the non-disturbed stationary state of the 
atom, we have 



^H(4))o = -^^(«)^L-;^> (26) 

where notation (...)o means that the expectations values are calculated using eigen- 
vectors of non-disturbed Hamiltonian Hq (p4[). If one neglects "deformation" of 
the atom by external gravitational field, then the electromagnetic field -^(a)(4) (see 



(|T9|), (14)) will affect only renormalization of the nucleus mass /i — > /i + A/i, where 



AyU = -j^- \J2^^)- Because of A/i//x <C 1 , we will not hereinafter consider this 

renormalization. 

Using perturbation theory, it is easy to show, that the required expectation value 
of the electromagnetic field tensor from equation (|20|) can be presented as Q: 

i^(a)(4) = (Fw(4)) = ^E [-'-^))+0{eUw), (27) 

where 0{e'^^,a-.) denote the terms quadratic in £[i)(j)- 

Direct calculation of quantity (J2[ ^^ ) ) is connected with mathematical 

difficulties in the framework of the perturbation theory. Therefore we shall use the 
following property of the hermitian operators. 

If observable L is described by the hermitian operator L, which does not depend 

o 

on time explicitly, the expectation value of the operator of evolution <L> of this 
observable L in the stationary state (for which the state vector is an eigenvector of 
the Hamiltonian) will be equal to zero: 

_<L>=<L>=0, (28) 



where 




cxpl 



1 r 



H,L 



(29) 



In general case the Hamilton operator of the atom in an external gravitational 
field always explicitly depends on proper time r. But, we can decompose the Hamil- 
tonian H of our quantum system in a vicinity of an arbitrary moment of time tq: 



H{t) = H{to) + 



m. 



(r - To) + ... 



'■TO 



and calculate the probabilities of transitions, caused by the term ( ^ ) (r — tq) 

in the operator H. If we will find, that during the time (r — tq) ~ At (where At 
is the characteristic atomic time) these probabilities are much less than 1, we will 
consider the eigenvector of the operator H as the quasistationary state of the atom. 
Numerical estimations show, that this quasistationary conditions for atom states are 
satisfied for a wide class of gravitational fields (such as the external gravitational 
field of the macroscopic black holes). In our case, the terms included in the Hamil- 
tonian Hi, such as 6 and ^(a), are very slowly change with time. Therefore, even 
in case of the ultrarelativistic motion of the atom, when the atom passes a short 
distances during a characteristic atomic time, the gravitational field can be consider 
as homogeneous with a high precision (except of the case of short-wave gravitational 
radiation). 

o 

We can introduce the operators of 3- velocity XA{ay 




H,X 



Aia) 



~ p 

m 

oo 



v(p) vM (0) v(t) 

A(a) - 3^^''(a){M)H(4)^A ^A ' ^{a){l3){T)^ ^A ' 



:CR(, 



and 3-acceleration XA{a)- 



(30) 



expl 



+ 



H, XA(a) 



1 I Ze XA(a) \ 1 



m \ p\ 



+ i; E 



m 






B^Bi^A 



2p 



\b 






{^A :Pa(0)} ,PA(a) 



+ 



+ 7::^^^^HP)ic.)ir)P^A^-W(^) 



m 



(31) 



By direct calculation it is possible to show, that, if the atom is in a quasista- 

o oo 

tionary state, the operators of 3-velocity XA{a) and 3-acceleration Xa(o) for each 
electron satisfy to the following equations: 



^lA(a)^=0, (xA(a)) = ^- 



(32) 



Taking into account the obvious explicit time- independence of the operators X 



(«) 



and -P4 and the fact, that they satisfy to following commutation relations: 



p(«) ^m 



^A '^B 



0, 



-S'-^^^'^Sab, 



p(a) p{/3) 
^A ^^B 



0, 



(33) 
(34) 



we receive from system of the equations ( P2D the following expression: 



^A{a) 



P\ 



-ZmW^""^- 



-(y2 [mc^A,ia) - ce('^)(^)(^)i?(,)(5)(4)M (S^ir) + LAir))] ) , (35) 

where we used, that "gravitational potentials" and ^(q,) very slowly vary with time: 

m(e(^)'j = %)(M)WH ({X«^),PH)}) - 0. (36) 

Because of we neglect of quadratic terms in £(i)(j), we can change from (...) to 
(...)q, where notation (•.•)o means that the expectation value is calculated in relation 
to eigenvectors of non-disturbed operator (24). Then we receive from the system 
( |l^ (15) the following expression for the equation of motion of the multielectron 
atom in an external gravitational field in the comoving reference frame: 



(// + Zm) W^(«) = -ce('")(^)(^)ii. 



mu)(a)ia) 



%)+ 



\ A I QJ 



"^C^ ( J2 ^A,{a) 



(37) 



\ A / 

The second term in the right part of the received equation (| 



-mc 



E' 



'A,(a) 



1 2 

—mc 

4 



R 



(4)(M)(4)(a)(r)- 



+%)(r)(4)(M)(a) + %){a)(4){T)(/.) J ( E ^/ ^A ) (^8) 

\ A / 

describes the interaction of an external gravitational field with the average quadru- 
ple moment of atom. From estimations made above, this term becomes actually 
small for atomic dimensions, if we assume, that the full spin of the atom (which 
consists of the spins of the nucleus and electrons, and also of the electronic orbital 
moment) is not equal to zero. The calculations of the expectation values of the op- 



erator J2 [ LA{a) + SA{a) I describing variation of the electron angular momentum, 
by means of (y) lead us to 



d_ 



^(a) + ( z2^^A{a) + SA{a)) 






\ A / 0-1 



(39) 



Finally, the equation of motion of the niultielectron atom it is possible to write 
down as: 



{^^ + Zm) Wm = -ce('^)(^)(-)i?MK.)(a)W ^(r) + 



'(r)- 



(40) 



+ { J2 i^Mr) + ^Air)) ) 
\ A /oJ 

In this way, we have demonstrated, using the some physically reasonable as- 
sumptions, that the motion of multielectron atom in an external gravitational field 
in a good approximation can be described by the system of classical Mathisson — 
Papapetrou equations 

(^ + Zm)^ = l-R^^.iu^e'^i^^UnJm, (41) 

^ = \u'Jn^, JmU"' = 0, (42) 

if we put as a classical angular momentum of atom the expectation value of the 
operator J^ of the total angular momentum of the system. In comoving rotating 
Fermi coordinates the quantity J has the following components: 

J{a) = S[a) + i 22(^A{a) + SA{a)) ) , J(i) = 0- (43) 

\ A / 

We note that if we also consider the influence of any external electromagnetic 
field y'^'F^^ on the atom motion, we have to add in the equation ( |4l] ) the following 
terms: 

^{Z + 1) (o)F^^n,- + ^ (o)f^ „;„£'""'=' Mfcwz, (44) 

where M^ is the 4- vector of the atoms full magnetic moment, which in comoving 
reference frame has the the following components: 

^^(") = ^(E(A.)+2^h)) +0^(")' ^A4)=0. (45) 

As it is visible from the received equation (|^), the trajectory of the atom will 
depend on the total angular momentum of the atom, which includes spins of the 
nucleus and electrons, and also orbital moment of electrons in atom. Hence, the 
motion of the atom in many respects depends on his quantum state. So, for example, 
the atoms, being in states with different mutual orientations of spins of the nuclear 
electrons, will move in general case along the different trajectories in an external 
gravitational field. Besides, the atoms, being in excited states with long life time, will 
be separated by an external gravitational field from atoms, which being in nonexcited 
states. Thus, the ultrarelativistic motion of atoms in the external gravitational fields 
can lead to a set of generally-relativistic effects, which can represent interest for 
astrophysics and allow us to check up the predictions of the general relativity on the 
quantum level. 
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